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Abstract. We prove weak-strong uniqueness in the class of admissible measure-valued solutions 
for the isentropic Euler equations in any space dimension and for the Savage-Hutter model of 
granular flows in one and two space dimensions. For the latter system, we also show the complete 
dissipation of momentum in finite time, thus rigorously justifying an assumption that has been made 
in the engineering and numerical literature. 
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1. Introduction 

A measure-valued solution to a paitial differential equation (or a system of equations) is, 
roughly speaking, a map that gives for every point in the domain a probability distribution of 
values, and that satisfies the equation only in an average sense. If this probability distribution 
reduces to a point mass almost everywhere in the domain, then the measure-valued solution is 
simply a solution in the sense of distributions. The main advantage of measure-valued solutions 
is the fact that, in many situations, they can easily be obtained from weakly convergent sequences 
of approximate solutions, even when the convergence of the approximating sequence to a distri¬ 
butional solution may fail due to effects of oscillation and concentration. 

Measure-valued solutions to hyperbolic conservation laws were introduced by DiPerna IIDiP85ll . 
He showed for scalar conservation laws in one space dimension that measure-valued solutions exist 
and are, under the assumption of entropy admissibility, in fact concentrated at one point, i.e. they 
can be identified with a distributional (entropy) solution. In other words, in this case the formation 
of fast oscillations, which corresponds to a measure with positive variance, can be excluded. 

In many other physically relevant systems, however, no such compactness arguments are avail¬ 
able, and existence of admissible weak (i.e. distributional) solutions seems hopeless. In such cases, 
the existence of measure-valued solutions is the best one can hope for. 

For the incompressible Euler equations, DiPerna and Majda IIDM87II showed the global exis¬ 
tence of measure-valued solutions for any initial data with finite energy. The main point of their 
work was to introduce so-called generalised Young measures, which take into account not only 
oscillations, but also concentrations. 

Subsequently, measure-valued solutions were shown to exist for further models of fluid and 
gas dynamics, e.g. compressible Euler and Navier-Stokes equafions IINeu931IKZ96ll or the Savage- 
Hutter avalanche model HGwiOSII . 

R G. is a coordinator of the International Ph.D. Projects Programme of the Foundation for Polish Science op¬ 
erated within the Innovative Economy Operational Programme 2007-2013 (Ph.D. Programme: Mathematical Meth¬ 
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2012/05/E/ST1/02218. The authors appreciate the support of the Warsaw Center of Mathematics and Computer 
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Measure-valued solutions have been criticised for being a too weak notion of solution. Indeed, 
the by now fairly standard procedure of establishing measure-valued solutions by viscous approx¬ 
imation, thereby circumventing delicate problems of compactness, suggests that the solution thus 
obtained does not carry enough information to be of much use. In particular, in the absence of 
admissibility criteria, measure-valued solutions are obviously non-unique to a large extent, as they 
only contain information on certain moments of the measure. 

It is therefore surprising that, in the case of the incompressible Euler equations, the so-called 
weak-strong uniqueness property was proved, on the whole space, for admissible measure-valued 
solutions by Brenier, De Lellis, and Szekelyhidi HBDLSl 11 . This means that if there exists a 
sufficiently regular (classical) solution, then every admissible measure-valued solution with the 
same initial data will coincide with the classical solution. Admissibility means that the kinetic 
energy of the solution never exceeds the initial energy. 

In fact, in BLio961 . P.-L. Lions required any reasonable concept of (very) weak solution to 
satisfy global existence and weak-strong uniqueness. For the incompressible Euler equations, 
therefore, admissible measure-valued solutions qualify. It is important though to emphasize the 
necessity of admissibility: Without this assumption, various examples are known where weak- 
strong uniqueness fails even for distributional solutions of incompressible Euler IISch93[[Shn971 
IDLS091IWiellll . Also, uniqueness need not hold for admissible solutions in the absence of a 
strong solution, see IIDLS101ISW121lDanl4ll . 

We consider in this aiticle two systems of equations in the realm of compressible fluid dynam¬ 
ics: The isenfropic Euler equations, 

dji -\- (^\\{hu) = 0 

dt{hu)-\-div{hu®u)-\-y{ kK^) =hG^ ^ ^ 

in any space dimension greafer or equal one, and fhe Savage-Huffer equations 

dt h -h div (hu) = 0 

0 ( 1 - 2 ) 
dt{hu) -|-div(/iM(8)n) + V(a/r) = h {—dB{u) -|-/), 

which make sense (from a modelling viewpoinf) in one or fwo space dimensions. Here, G and 
/ are exfernal force densifies, and B{u) is a maximal monotone sef-valued map. The lafter sys- 
fem describes fhe evolufion of fhe depfh-averaged velocify and heighf of some maferial sliding 
over an inclined slope. The maferial is subjecf fo fhe so-called Coulomb-Mohr friction law. For 
comprehensive sfudies, including derivafion, numerical compufafions and experimenfal resulfs on 
sysfem (11.21) and ifs various modificafions, we refer fo IISH89IIGWH99BBW04IIBMCPV03IIGTN031 
IGC07[|HWP051lPBM081IZPTN10l . among ofhers. We prove (cf. Theorems l4.2l and l5.2l belowl: 

Theorem 1.1. Let (H,U) be a solution of dEB or (fOI) such that H is Lipschitz continuous in 
[0, T] X T” and U € C^([0,r] x T”). Assume also H > c > 0 for some constant c. Then every 
admissible measure-valued solution of (EB or (fOl) . respectively, with the same initial data as 
{H,U) coincides with {H,U). 

Of course, fhe precise definitions of measure-valued solufions and admissibilify will be given 
below. If should be menfioned fhaf weak-sfrong uniqueness for admissible measure-valued solu¬ 
fions was proved in IIBDLSllll for general hyperbolic systems of conservation laws, but this was 
done only for oscillation measures. 

Moreover, the results in IIBDLSllll are valid even for Lipschitz continuous strong solutions. 
Owing to commutator estimates analogous to the ones provided by Constantin, E and Titi for the 
incompressible Euler system in IICET94II . the result of Theorem I 1.1 l ean be obtained assuming only 
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Lipschitz continuity of U, and Sobolev regularity of H. In fact, as in HBDLSllll it is sufficient to 
assume only that the symmetric part of Vf7 be bounded. We omit details. 

Weak-strong uniqueness for compressible Euler models appears important in the light of several 
recent examples of non-uniqueness of admissible weak solutions HDLS 101IChil41IOC141ICFK151 
ICDLK14llFeiT4l . For the Savage-Hutter equations, such examples were very recently constructed 
in IIF(1S(115II . 

To prove the weak-strong uniqueness, we follow the general strategy of IIFJN12II (where weak- 
strong uniqueness was proved for the compressible Navier-Stokes equations), but we have to adjust 
these arguments to the measure-valued framework. For the Savage-Hutter system, an additional 
issue is to give a meaningful definition of measure-valued solutions that accounts for the multi¬ 
valued nature of the fricition term B{u) in (11.21) . Such a definition was proposed in lIGwiPSH and 
we will use it here as well. 

If the force / is time-independent and ||/||oo < d, then a special class of solutions to (11.21) is 
given by M = 0 and h independent of time and such that 


Wh{x) - 


fix) 


2a 


< — for almost every x. 
2a 


Observe that our weak-strong uniqueness result allows to take {H,U) as such a stationary solution, 
so that in particular every such solution enjoys uniqueness in the class of admissible measure¬ 
valued solutions. 

For the Savage-Hutter model we also prove the following result: 


Theorem 1.2. There exists a finite time Q <T < oo, only depending on the parameters in (11.21) 
and the initial data, such that every admissible measure-valued solution of (fO) starting from 
such initial data has zero momentum for almost every time t > T. 


In particular, this implies that every admissible weak (distributional) solution becomes station¬ 
ary after finite time. This highlights the importance of stationary solutions as well as the role 
played by the admissibility condition: Indeed, in IIFGSG15II non-admissible weak solutions were 
constructed whose momentum does not decrease to zero. The result is a rigorous justification of 
empirical and numerical observations of deposition of material after finite time, IISH91llFNBB+08l 
ICGM12I . The finite-time runout of solutions is essentially used at the modelling stage as providing 
data for calibration of the system. This property was assumed in numerical simulations, however, 
to our knowledge, never proved. 

Fet us remark that for the one-dimensional Savage-Hutter model, we obtain a fairly complete 
picture: Existence of admissible global in time weak solutions is known IIGwi02ll . they enjoy 
weak-strong uniqueness, and become stationary after finite time. 

Similarly, for the compressible Euler system in the one-dimensional case there exist global in 
time admissible weak solutions having the weak-strong uniqueness property, see IIDiP831IFPS961 
IFPT94I . 

Finally, let us point out some difficulties in extending our results to other domains than the 
torus. On the whole space, we can no longer require the denstity H to be uniformly bounded away 
from zero and the initial energy 

f I 2 2 

/ -hofiol -\-ahQdx 
JR" 2 

(and similarly for (11.11) 1 to be finite at the same time. On domains with physical boundaries, 
however, we do not even expect weak-strong uniqueness to hold, since a counterexample has been 
exhibited in IIBSW14II in the case of the incompressible Euler equations. 
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2. Notation 

We fix here some notation that will be used throughout the paper. 

The n-dimensional torus will be denoted by T" := 

Let n C M” be a measurable subset or H = T". The set of locally finite nonnegative measures 
on will be denoted If X is a measurable subset of M"*, then will be the set of 

probability measures on X. 

Let m G The space ^{X)) is then defined as the space of maps v : —)• 

I^{X), X i-A Vx, which are weakly* measurable with respect to m; that is, for every (p G Cc{X) the 
map 

Vi-A / q){X)dVx{X) 

Jo. 

is m-measurable. If m is Lebesgue measure, we simply write ,^(X)). If Q. has the form 

[o,r] X n for some measurable subset C M", then dx denotes n-dimensional Lebesgue measure 
and dt one-dimensional Lebesgue measure. The Dirac mass centred at x will be denoted as 8x, as 
usual. 

The ni-dimensional unit sphere is written With D we mean the topological closure of a 
subset of M”. We write M+ for the set of non-negative real numbers. 

In the case D = [0, T] x fl, we will use measures of the form m = nij (g) dt\ this means that, for 
every set of the form z xU, where T C [0, T] and 17 C D are measurable subsets, 

m{zxU)= J mt{U)dt. 

The differential operators V and div are applied only to the spatial variables. If u and v are 
vectors, then n (g) v denotes the matrix defined by {u (g) v)/y = UiVj. The divergence of a matrix field 
is understood to be taken row-wise. 

Further notation will be introduced as we proceed. 


3. Generalised Young Measures 


We recall briefly the notion of generalised Young measures, which were introduced by DiPerna 
and Majda IIDM87II and refined by Alibert and Bouchitte IIAB97II . Further details can be found 
e.g. in BKRlOllSWnil . 

Young measures are used to represent weak limits of nonlinear functions of weakly convergent 
sequences. More precisely, suppose Q. C ffi" is a measurable set or D = T" {n > 1), and is 

a sequence of maps bounded in {m > 1). 

Then it was proved in IIAB97II that there exists a subsequence (not relabeled) as well as a 
parametrised probability measure v G iW")) (which is identical with the ’’classical” 

Young measure), a non-negative measure m G and a parametrized probability measure 

v“ GL“(a,m;,^(S"'-i)) such that 

f{x,u„{x))dx ^ [ f{x,X)dVx{X)dx+ ( /“(x,j3)dv”(j8)m 

Jkp' js™-! 


weakly* in the sense of measures. Here, /: D x —)■ M is any Caratheodory function (measurable 

in the first and continuous in the second argument) whose recession function 


/“(x,j3) := lim 

/S'^/3 


s 


is a well-defined and continuous function on D x S™ ^. Note that such an / will have at most 
linear growth. If its growth is sublinear, then /“ = 0. 
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Notice also that v," is only defined m-almost everywhere. 

In order to properly define measure-valued solutions to compressible fluid equations within the 
framework of Alibert-Bouchitte, we need a slight refinement which allows us to treat sequences 
whose components have different growth. Let {uk,Wk)k be a sequence such that {ut) is bounded 
in LP{Q.-,MJ) and {wk) is bounded in (1 < p,q < oo). Define the “nonhomogeneous unit 

sphere” 


■= {(Pulii) G + |/32|2‘? = 1}. 

Then, there exists a a subsequence (not relabeled) and measures v G L“(D; m G 

v“ G such that 

f{x,u„{x),Wn{x))dx^ [ f{x,Xi,X2)dVx{XuX2)dx+ [ /“(A,j3i,j82)r/v“(j8i,j82)m 

JK'+m 


in the sense of measures; this is valid for all integrands / whose /^-(^-recession function exists and 
is continuous on D x The /7-<7-recession function is defined as 


/“(x,j3i,j32) := lim 


fix',s‘ip[,sPI5') 

sP‘1 


The case p = 2, q = I was treated in Subsection 2.4.1 of IISW12II and the extension to general p 
and q is straightforward. 

Let us quote another fact which is important for measure-valued solutions of time-dependent 
equations with bounded energy: If D = [0,r] x D for some measurable D C M" (or D = T"), 
and if the sequence (m„,w„)„ is bounded in L°°([0,r];LP(D) xL'?(D)), then the corresponding 
concentration measure m admits a disintegration of the form 


m = mt{dx)®dt, 

where 1 1 —m, is bounded and measurable viewed as a map from [0, T] into The proof of 

this statement was given in IIBDLSllll . 


4. Weak-strong uniqueness for measure-valued solutions of the 

COMPRESSIBLE EULER EQUATIONS 
We consider the compressible Euler system 
dth + di\(hu) = 0 

(4.1) 

dt{hu)+(^i\{hu®u)-\-'W{Kh^) =hG. 

Here, /i: [0,r] x T" —>■ M, n : [0,r] x T" —)• M", and G : [0,r] x T" —)• M”, /> 1. We set the constant 
K > 0 equal to one in order to save some writing, remarking however that all computations remain 
unchanged for general K. 

4.1. Measure-valued solutions. We apply the abstract framework from the previous section, with 
I = I, m = n, p = j, and q = 2, in order to define the notion of measure-valued solution of (14.11) . 
Consider a generalised Young measure 

(v,,„m, ) G ([0, T] X r; ^(M+ x M")) x ^+{[0, T] x T") x ([0, T] x T",m; ^(§+)) , 

where we wrote 
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We will use the variables {X\,X') £ x M" and (j3i,j3') £ as dummy variables when inte¬ 
grating with respect to Vt^x and , respectively. One should think of Ai,j8i as representing h and 
X',p' as representing Vhu. We also use the common notation 

(F(A,,A'),Vm) := [ F(A,,A')r/v,,(Ai,A') 

v/R+xR" 

and analogously for v”. 

If we consider a function /: [0, T] x T” x M+ x M" —)• M which has an appropriate 7-2-recession 
function as defined in Section |3l we use the shorthand notation 


f{dtdx) := {f{t,x,-,-),v,^)dtdx+{f°{t,x,-,-),v",,)m{dtdx). 
For instance, we have 

^ = (Ai,v) 

hu = , v) 

hw^u = {X' ® X', v) + (8) j3^, v“)m 

W=(|Af,v) + (|j8f,v>i 
(AiG,v) =^G. 


We say that (v,m, v“) is a measure-valued solution of (14.11) with initial data {ho,uo) if for every 
T£ [0,r], i//£C‘([0,r] XT";]R), (p£Ci([0,r] xT';M") it holds that 


0 JT" 


dt\lfh + V\(f ■ hudxdl + / i/r(A,0)/io — >/r(A, T)(fA = 0, 


/7 

Jo Jt" 


/T" 


dt<J) ■ hu + W(p \ hu® u + div (ph^ — (p ■ hGdxdt 


+ / (p{x,0) -hoUQ — (p{x,T) ■hu{x,T)dx = 0. 
JT" 


(4.2) 


It is part of the definition that all the integrals have to exist for any choice of test functions, in 
particular for the initial data we require /iq £ hoUQ £ L^. 

Let us set 


1 - 


1 


Emvsit) ■= / -h\u\^{t,x)-\ - -hy{t,x)dx 

JT" 2 7 — 1 

for almost every t, and 


En- = 


f 1 1 

/j,„ 2^o\uoHx) + Y^hl{x)d. 


We then say that a measure-valued solution is admissible if 

Emvs{t)^EQ+ [ [ hG ■u{s,x)dxds (4.3) 

JO JT" 

in the sense of distributions. An elementary computation yields the well-known fact that the energy 
is conserved for smooth solutions (i.e. (14.31) holds with equality), whereas the inequality becomes 
strict upon the formation of shocks. 


Remark 4.1. The global existence of measure-valued solutions for (14.11) was proved by Neustupa 
in IINeu93ll . However he used a different formulation of the Young measure, as the formalism of 
Alibert-Bouchitte BAB9711 was not yet available. One can however rewrite the solutions of BNeu93B 
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in the form presented here. Neustupa’s solutions can be seen to be admissible, as they can be 
obtained e.g. from an aiTificial viscosity approximation. 

4.2. Weak-Strong Uniqueness. 

Theorem 4.2. Let G G L“([0,r];L2(T”)) and suppose H G Wi’“([0,r] x T"),U G C'([0,r] x 
T”) is a solution of (14.11) with initial data hn >c>0,ho gL^(T”), /io|moP GL'(T”), andH{xd) > 
c > 0for some constant c and all {t,x) G [0, T] x T”. If^vpn, v") is an admissible measure-valued 
solution with the same initial data, then 

= 0 - 

Proof. Let us first define for a.e. t G [0,r] the relative energy between {H,U) and the measure¬ 
valued solution as 


Erei(t) = [ lh\u-U\^ + -^hr - —Hr-'^h + mdx 

Jt' 2 7-1 7-1 

pi 2 1 JT' 

7t" 7- 1 7- 1 

1 


/ 

JT' 


+ 


/T" 7" 




Here, the measure nit G (TT”) is obtained by the disintegration m{dtdx) = {dx)®dt, which is 
well-defined thanks to the admissibility (cf. Section[3l). Note that the strict convexity of | • |^ implies 
that the relative energy is always non-negative. Then it is straightforward to observe that Erei (t) = 
0 for a.e. t implies Theorem l4.2l Indeed, defining the projection operators : (Ai, A') i—)■ Ai and 
K^' : (Ai,A') I—)• A', we observe that the strict convexity of | • |^ implies that Ti^^Vt^x = for 
a.e. t,x and hence Vt^x = S(^H{t,x)) ® ^^''^t,x- Using the first term in the relative energy allows to 
conclude 

Setting (p = U in the momentum equation (the second equation of (14.21) 1. we obtain 


hu-U{x)dx= ho\uo\^dx-\- 


+ 


to JT” 

WdivU + hG-Udxdt. 


hu ■ dfU -\-hu®u\ VU dxdt 


(4.4) 


0 JT” 


Similarly, setting i/r = j\U\'^ and then i/r = yW ^ in (14.21) yields 
1 


2 JT" 


\U{t)\ h{T,x)dx = 


0 JT" 


U ■ dfUli VUU ■ hudxdt 


Ik 


uoph^dx 


(4.5) 


and 


/T" 


yEl^ ^{z)h{'c)dx = 


0 JT" 


7 ( 7 — 1 )//^ ^dtHh-\-y(y—\)H'^ H■ hudxdt-\- / yh^dx, 


T" 


(4.6) 


respectively. 

Next, we can write the relative energy as 


Erel{T:)= f 

JT" 2 7 — 1 


1 


h^dx-y I H^dx-\-— I \U\^hdx — 

T" 2 JT" 


U■hudx— 


T" 


IT” 7 — 1 


-H^-^hdx 


1 


= Emvs{'^)+ I H^dx + - I \U\'^hdx— I U-hudx — 

/T" 2 ” 


/- 

JT" 


T" 7 — 1 


-H^-^hdx 
















PIOTR GWIAZDA, AGNIESZKA SWIERCZEWSKA-GWIAZDA, AND EMIL WIEDEMANN 


(all integrands evaluated at time t). Next, using the balances (14.41) . (14.51) . (14.61) for the last three 
integrals, we obtain 


EreliT:) = + [ W^dx 

JT' 


+ [ [ U ■ dtUh + VUU -hudxdt + f -\uo\^h(jdx 
Jo JT" JT" 2 

— / ho\uo\^dx— / / hu-dtU+hu®u'.WUdxdt 

JT" Jo JT" 

— [ [ div U — hG ■ Udxdt 
Jo JT" 

■ hu)dxdt — 


)o JT" 

and using (15.41) we have, for a.e. T, 

Erei{'^)^— [ hl.dx+ [ [ hG-u+ [ H^dx 
Jt" Jo JT" JT" 

+ [ f U ■ dtUh + WU ■ hudxdt 
Jo JT" 

rt r _ _ 

hu ■ dtU + hu®u'. V[/ dxdt 


V 7 — 1 


^ hj^dx, 


to JT" 
r'c 


hrdivU -hG-Udxdt 


to JT" 

(•T 


to JT" 

Next, we collect some terms and write 

J H^dx— j hj^dx — 


( 7 //^ ^dtElh + yEl^ ^VH-hu)dxdt. 


T" 

^ d 


to JT" 
y-2: 


yH^ ^dfHhdxdt 

—H^-yH^-^d,Hhdxdt 
dt 

r [ yH ^-' dtH - yH^-^dtHhdxdt 
Jo JT" 

yH'^-^dtHiH -h)dxdt, 


to JT" dt 

rt 


rt 
to JT" 


(4.7) 


(4.8) 


and 


hG ■ u — hG ■ Udxdt 


0 JT" 


[ [ hG ■ {u — U)dxdt, 
Jo JT" 


0 JT" 


U ■ dfUh + yUU ■ hu — hu ■ dfU — hu<S> u : VUdxdt 

[ [ dtU ■ h{U — u) d-yu : hu^ {U — u)dxdt. 
Jo JT" 


(4.9) 


(4.10) 


Indeed, the last two equalities can be verified by writing the expressions in the ’’coarse-grained” 
overline notation explicitly in terms of the Young measure (v,m, v"). 
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Plugging equalities (14.81) . (14.91) . (14.101) into (14.71) . we arrive at 

Erei{T:) < r [ -h)dxdt 

Jo Jt" 

- [ [ hG ■ {u — U)dxdt 
Jo 7 t " 

[ [ dtU ■ h{U — u)+ ^11hu® {U — u)dxdt 
Jo Jv 

— [ [ hy diiwUdxdt — [ [ ■ hudxdt. 

Jo Jt' Jo Jt" 

For the last two integrals, we have by the divergence theorem 

— f f hy divUdxdt — f f ■ hudxdt 

Jo Jt" Jo Jt” 

pr r _ _ 

-hydiv U + yH^ ^VH ■ {HU - hu) - yH^ ^VH ■ HUdxdt 


10 JT" 


0 JT" 


{H^ - hy) div U + yH^ ■ {HU - hu)dxdt. 


Inserting this back into (14.111) and observing that, by the mass equation for {H,U), 

yH^-^dtH{H - h) + yH'^-'^divUH{H -h) + yH^ \H • HU = yH^-^U ■ VHh, 

we get 

Erei{T:) < r [ ■ VHh{U - u)dxdt 

Jo Jt' 

+ [ [ hG-{u — U)dxdt 
Jo Jt” 

+ [ [ dtU ■ h{U — u) d-VUhu® {U — u)dxdt 
Jo 7 t " 

- r [ yH^^^divU{H-h)dxdt+ f [ {H^-hy)divU. 

Jo Jt" Jo Jt‘ 

The expression in the third line can be rewritten pointwise as 


dtU ■ h{U — u) + yU : hu®{U — u) 


dtU-h{U-u) + VU : U®h{U-u) 


= —{dt{HU)+div{HU®U))-h{U-u) 


H 


= G-h{U -u)- yH^ ^VH ■ h{U - u). 

Putting together (14.131) . (14.151) . and (14.161) . we obtain 

Erel{^)<C\\U\\cl fErel{t)dt. 

Jo 

Finally, from Gronwall’s inequality it follows that Erei{'c) = 0 for a.e. t. 


(4.11) 


(4.12) 


(4.13) 


= dtU ■ h{U -u) + yu -.hU ®{U-u) + yu ■. h{u -U)®{U-u), 
and the integral of the last term as well as the last line in (14.131) can both be estimated by 

C||f/||c. fErel{t)dt. 

Jo 

For the remaining terms in (14.141) we obtain, using the momentum equation for {H,U), 


(4.14) 


(4.15) 


(4.16) 


□ 
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5. Savage-Hutter system 
We consider the two-dimensional Savage-Hutter model 
dth + d\\r{hu) = 0 

o (5.1) 

dt{hu) +d.\y{hu®u) + W{ah^) = h {—dB{u) + f). 

The one-dimensional case can be treated similarly. Here, /i: [0, T] x —)• R, n : [0, T] x —)• R^, 
/ : [0, r] X —)■ R^, and a > 0 and d > 0 are constant. By B{u) we denote the subdifferential of 
u I—)■ \u\, so that B{u) is multi-valued such that 

»(«)=(!!_ 

[Bi(0 ) ifM = 0. 

Consequently, the equality sign in the second line of (15.11) should really be an inclusion. We will 
stick however to the formulation (15.11) . thereby slightly abusing notation. 


5.1. Stationary solutions. If in (15.11) / is independent of time, then a special class of solutions is 
given by M = 0 and any h = h{x) > c such that 

V/i —— < — for a.e. x. 

2a 2a 


5.2. Measure-Valued Solutions. We recall the notion of measure-valued solution of (15.11) from 
IIGwi05ll in the notation used therein (in fact, there the problem was treated on the whole space, 
but it can easily be adapted to the torus). The author considers system (15.11) with a right-hand side 
given by hf{x, y/hu), where 

/(a, Vhu) = -d -h / 

and for u = 0 the mapping / takes values in the closed unit ball. To handle this multi-valued 
(monotone) term, let us first recall from IIGwi05ll the following observation, see also IIGZG071 
BGMSGf)^ for a similar approach. 

Lemma 5.1. Let / : R" —R" (/ : R” —)• 2*") be a monotone function (monotone mapping). Then 

(f + ld)-^ :R"^R" 


and 

/o(/ + W)-i :R'’^R” 

are Lipschitz functions. Above we understand /o (/ + Id)^^ = Id — Ido (f Id)~^ 
Moreover, for any continuous function g : R” ^ R”, 


go(/-yW)-' :R”^R" 


is a continuous function. 

Under the assumptions that (/io,/iomo) C {^\uo\^ho + aQig)^'^ dx < o°, there exists 

a triple of measures 

eL“ ([0,r] X t 2;'^(R+ x R^)) x .^+([0, T] x T^) x L;; ([0,r] X T^m;'^(S+)) , 
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such that 


/ / hdt(pi+m-'Vx(Pidxdt = / hQ(pi{0)dx, 

J[0,7')JK2 JR2 

mdt(p2dx+(e+ dp,'Vj((p2) — / f(p2dx>dt= / hQUQ(p2iO)dx 
l[0,T) Lv/R2 JR2 J JR2 

for all (pi,(p 2 £ C'c([0)7’) x and for almost all t G [0,r) it holds that 

{Tr {e {t)) + dp - [ \l:\uo\^ho + a{ho)^\dx< [ xdxdt, 
dR2 [2 J J[o,r)xR2 

where 


(5.2) 


h{t,x) = Ai dpt,xi^), 

2R_|_ xR2 

p(t,A) = f Af dpt,xi^) + ( / Pi r/v“ 

m(t,A) = [ x/xli-f + Idy\x,{X2,X3)) dpt^xi^), 

JR+XR2 

e(t,A) = [ {-f + Idy\x,{X2,X3))0{-f + Idy\x,{X2,X3))dpt,x{^) 

JR+xR2 

+ (^2 (/32,^3) <8* (/32,^3) dv“ 

Tr(e)(t,3c) = [ {-f + Idy\x,{X2,X2)) ■ {-f + Idy\x,{X2,X3)) dptA^) 

JR+XR2 

+ (^2 ^2 + ^3 dVfAP)^>n, 

f{t,x) = f Xj o {-f + IdA\x,{X2,X3))dpt.x{X), 

2R+xR2 

x{t,x) = [ Xif O {-f + Idy\x,{X2,X3)) ■ {-f + IdA\x,{X2,X3))dptA^) 

JR+XR2 

for almost all (t,x) G [0, T] x T^. 

Define for every /r-measurable set A x B C x the push-forward of p through the map 
{-f + Id) as 

v(A xB) := {-f + Id)#p{AxB)=p{A x {-f + IdA\B)). 

Hence using again the variables {X\ ,X') G x (which correspond to X\ , A 2 , A 3 ) and (j 8 i, j 8 ') G 
(corresponding to pi,p 2 ,P 3 ), when integrating with respect to v,_i and v,", respectively, the 
problem can be translated to 

/i= (Ai,v) 

lfi = {Xl,v) + {PlvAm 

hu = {x/xlX' ,v) 

hw^u = {X' ® X', v) + {p' ( 8 ) p', v“)m 
lAA = {\xf,v) + {\pf,vAm 
h{-dB{u)+f ) = (Ai/o {-f + Id)-\x,X'),p). 
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We say that (jU,m,/i") is a measure-valued solution of (15.11) with initial data (/jo,mo) if for every 
TG [0,r], i//GC‘([ 0 ,r] xT 2;M), (pGCi([ 0 ,r] xT2;R2) it holds that 


/o JT 2 

rT 


dt^fh + V^f-hudxdt + / \ff{x,0)hQ — \ff(x,T)h{x,'c)dx = 0, 
JT 2 


0 JT 2 


dt^ ■ hu + Vq): hu<Siu + adiv (ph^ + ip • h{—dB{u) + f)dxdt 


+ 


/ (p{x,0) ■ Iiquo — (p{x,T) ■ hu{x, 
Jt 2 


t) = 0. 


For a.e. t, we set 


Emvsit) ■= / -h\u\^{t,x)+ah^{l,x)dx 
Jt- 2 


and 


Eo:= [ ^ho\uo\^{x)+hl{x)dx. 

JT" 2 

We say that a measure-valued solution is admissible if 

Emvs{t)<Eo-f [ dh{B{u)-f)-u{t,x) 

Jo JT" 


(5.3) 


(5.4) 


in the sense of distributions. 


5.3. Weak-Strong Uniqueness. 

Theorem 5.2. Let / G L“ ([0, T]; (T^)) and suppose // G W ‘ ’“ ([0, T] x T^), U G C' ([0, T] x T^) 
is a solution of (15.11) with initial data /iq > c > 0, /iq G L^(T^), /jo|moP S L' (T^) andH{t,x) > c > 0 
for some constant c and every (t,x) G [0,r] x T”. If {v,m,v°°) is an admissible measure-valued 
solution with the same initial data, then 


^ = 0 - 


Proof Let us first define for a.e. f G [0,7] the relative energy between {H,U) and the measure¬ 
valued solution as 


Erel{l) = [ -h\u — U\^-\-a{h—H)'^dx 
jt 2 2 

= / ,Vt^)dx+l- [ {\l5'\^,v^f/dmt{x) 

Jt 2 2 2 Jt2 ' 

+ a [ {\Xi-H\'^,Vt^x)dx-\-a [ {jif,vff)dmt{x). 

JT2 ,Jj2 


Here, the measure G ./#+(T^) is obtained by the disintegration m{dtdx) = mt{dx)^dt, which 
is well-defined thanks to the admissibility. 

It is straightforward to observe that Erei{t) = 0 for a.e. t implies Theorem 15.21 
Following the computations of Section |4] we arrive at 


Erei{^)+ r [ h(dB(u)-dB{U))-(u-U)dxdt <C\\U\\ci f Erei{t)dt. 
Jo Jtp Jo 


Finally, since B is monotone, the integral on the left hand side is non-negative, and from Gronwall’s 
inequality it follows that £'re/(T) = 0 for a.e. t. □ 
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6. Dissipation of Momentum in Finite Time 

Theorem 6.1. Let (v,wi,v") be an admissible measure-valued solution of the Savage-Hutter 
system (EB with initial energy Eq and 

II/IIl“(r+xT2 ) < d- 

Then there exists 0 <T <oo such that 

M{t) := / h\u\{t,x)dx = 0 for almost every t > T. 

Moreover, there exists a constant C depending only on d — ||/||oo and a such that 
T < 

Proof For the momentum we have 

M{t)= [ h\u\{t,x)dx= [ {\/Xi\X'\,Vtx)dx. 

Note in particular that the momentum does not concentrate, i.e. the 2-2-recession function of 
(Ai,A') 1 -^ Ai|A'| is zero. 

For the following estimate, we use Jensen’s inequality applied to the function | • (recall that, 
according to our convention, the torus has measure 1), then Young’s inequality, 

, aP b‘1 \ \ , 

ab <-1-, —I— — 1, 

p q p q 

with the conjugate exponents 3 and 3/2, and finally the admissibility assumption: 

/ /■ \ 

-t/ i^i^^t,x)dx + ^ f {\X'\^,Vt,x)dx 

J Jt^ .5 Jfi 

<C{a)E^,,{t) <C{a) (^Eo-- \\f\\^){s/I'i\Xlv,,x)dxds^ 

= C{a) (d- ll/lloo)^ M{s)ds^ , 

where 

C(»)=max|l 1/ 

Therefore, for almost every t, M{t) is less than or equal to the solution of the integral equation 

= C{a)EQ - C{a){d - ||/||„.) f M{s)ds 

Jo 

or equivalently (after differentiating) 

Mft) = -^C(a)(r/- ||/||„)M(^)"/^ M(0) = {C{a)Eof^\ 

The solution of this ordinary differential equation is easily computed as 

M^t) = I - IC{a){d - ||/||.)f)]' if 3{C{a)Eoy/^ - |C(a)(d - ||/||.)t > 0, 

10 otherwise. 
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In fact, [5(3(C(a)£'o)*/'* — ^C{a){d — ||/||<»)t)]^ would also be a solution for all times, but we 
know a priori that M{t) must be non-negative. 

It follows that there is a time 

after which M{t) = 0 almost everywhere. 

□ 


Corollary 6.2. Let [h, u) be an admissible weak solution of the Savage-Hutter equations with 
initial energy Eq and ||/||oo < d. Then there exists a time Q <T < oo such that for almost every 
t >T, (/r, u) is stationary, i.e. u{t,x) = Ofor almost every t >T and x £ T^, dth{t,x) = 0/or almost 
every t > T, x (z T^, and 




Moreover, T satisfies the estimate ofTheorem \6.1\ 


Proof As every admissible weak solution can be viewed as an admissible measure-valued solution 
via the identification v = 5^^ m = 0, from Theorem ib.ll we obtain a time T such that after this 
time, the momentum is zero: 

/ h\u\dx = 0 for almost every t > T. 

7t2 

Therefore, the Savage-Hutter equations reduce to dfi = 0 and 
V{ah^) = —dhB{u) + hf. 

The latter is clearly equivalent to \Wh — f /2a\ < d/2a, given that \B{u)\ <1. □ 
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